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Abstract 

During thunderstorms on 2008 September 20, a simultaneous detection of 7 rays and electrons 
was made at a mountain observatory in Japan located 2770 m above sea level. Both emissions, 
lasting 90 seconds, were associated with thunderclouds rather than lightning. The photon spectrum, 
extending to 10 MeV, can be interpreted as consisting of bremsstrahlung 7 rays arriving from a 
source which is 60 — 130 m in distance at 90% confidence level. The observed electrons are likely 
to be dominated by a primary population escaping from an acceleration region in the clouds. 

PACS numbers: 82.33.Xj, 92.60.Pw, 93.85.-q 
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INTRODUCTION 



Energetic radiation bursts from thunderstorm activity have been observed by space obser- 
vatories as well as ground-based detectors. Interestingly, these bursts appear to be classified 
into two types according to their duration. One comprises short bursts lasting for milli- 
seconds or less, mainfy^bserved from the upper atmosphere of the Earth O], natural 



lightning discharges 



mospner 



Though not necessary homo- 



^, and rocket-triggered ones 
geneous, they generally occur in association with lightning discharges. The other consists 
of long bursts with duration of a few seconds _to_ a few minutes, mostly detected at aircraft 
at high mountains 1^, 



altitudes [8| 
Sea 15 
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I2I . [13I . Q , and around the coastal area of Japan 



16|. Unlike the short bursts, few of them appear to be clearly associated with 
lightning discharges. 

To date, at least some of these radiation bursts are thought to result from relativistic 
electrons, which are produced in turn by a mechanism involving runaway and avalanches of 
electrons in the atmosphere. The mechanism, first proposed by Gurevich et al. ITj, needs 
intense electric fields, high-energy seed electrons, and spatial length long enough for the 
avalanches to develop. When these conditions are fulfilled, the electrons can gain energy 
from the electric fields fast enough to exceed their ionization loss. Then, they will be 
accelerated to relativistic energies, and produce bremsstrahlung 7 rays. 

Actually, 7-ray spectra extending to ~10 MeV have been observed from both short and 



long bursts 



IGj, with their properties generally consistent with the bremsstrahlung 



emission scenario. However, according to Monte Carlo simulations considering a feedback 
mechanism in which back-scattered photons and/or positrons successively supply seed elec- 
trons in the high-electric-field region, electron avalanches are likely to last only for microsec- 



onds 



3, Q, 



20|]. Observationally, only a few examples of long bursts have so far been 



reported to extend to MeV energies. Therefore, it is not yet clear whether the runaway 
electron avalanche scenario, which is generally successful on some short bursts 

0.9 

, can 

also explain the prolonged events. Furthermore, we need to explain the fact that some 

y, Q, Q, Q, H, Q, while others involve only 



long bursts consist so^ 
charged particles [13 



fily of x/7 rays 

14| . Here, utilizing a detector that is capable of separating photons 
from charged particles, we first report on a simultaneous detection of 7 rays and electrons 
in one long burst, with the photon spectrum extending to 10 MeV. 
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EXPERIMENT 



The Norikura cosmic-ray observatory (36°6'N, 137°33'E), belonging to the Institute for 
Cosmic Ray Research, the University of Tokyo, is located on a high mountain in Gifu 
prefecture, Japan. Due to its high altitude, 2770 m above sea level, as well as meteorological 
conditions from summer to autumn, the observatory is frequently shrouded in thunderclouds. 

Outside the observatory, wc installed a radiation detection system on 2008 September 4, 
and operated it until 2008 October 2. The system consists of a spherical Nal scintillator 
with a diameter of 7.62 cm, and a 45 x 40 x 0.5 cm^ plastic scintillator which is enclosed in 
an aluminum box with the top and bottom being 1 mm and 3 mm thick, respectively. Each 
scintillator has a photomultipher (HAMAMATSU R878) of its own, and each output signal is 
fed to a self-triggering electronics system, incorporating a 12 bit analog-to-digital converter 
(ADC). The detected events are collected in the order of arrival without individual arrival- 
time information, and recorded every 1 sec synchronized by Global Positioning System. 
Thus, both scintillators have a time resolution of 1 sec. The Nal scintillator was operated 
over 10 keV —12 MeV, while the plastic scintillator, placed above it, was operated with a 
threshold energy deposit of 500 keV. 

The plastic scintillator has a high sensitivity for charged particles while it is much less 
sensitive to photons. Therefore, when the Nal detector is triggered by an incident particle 
without a coincident (within 0.7 /^sec) signal in the plastic scintillator, the incident particle 
would be a photon, not a charged particle. Thus, the plastic signals mainly represent 
incoming charged particles including electrons, and using them in anti-coincidence, the NaT 
signals can be considered to indicate 7 rays with little (~ 0.1%) contamination by electrons. 

Aiming at an independent radiation observation, we installed, inside the observatory, an- 
other inorganic scintillation detector. This indoor detector uses a cubic bismuth germanate 
(BGO) scintillator of 5 x 5 x 15 cm^, which is coupled to a photomultipher (HAMAMATSU 
R7600U-203), and observes in the 100 keV — 5 MeV range. The signal is pulse-height an- 
alyzed by a self-triggering 8 bit ADC, and the pulse height is recorded, on event-by-event 
basis, with arrival-time information of 20 ms resolution. 

In addition to the radiation observations, we measure optical and electric field variations 
outside the observatory. The visible-light intensity is measured by an optical sensor using a 
photodiode (HAMAMATSU S1226-8BK) having a sensitivity to photons with wavelength of 



3 



320 — 1000 nm, with its peak at 750 nm. Output signals of this optical sensor are fed to a 12 
bit ADC every 1 sec, and recorded as voltages between and +10 V with a resolution of 5 
mV. The electric field is measured by a commercial electric field mill (BOLTEK EFM-100). 
Its outputs are also collected by a 12 bit ADC every 1 sec, and converted to the electric field 
strength between ± 100 kVm~^ with a resolution of 50 Vm^^. 

RESULTS 

Examining the data over a period of September 4 to October 2, we found a long-duration 
enhancement during thunderstorms on September 20, but found no other events of similar 
properties. Figure [1] shows count histories of > 10 keV NaT and > 500 keV plastic scintil- 
lators, between 15:15 and 16:15 UT on September 20. Besides gradual count-rate changes 
mainly due to Radon rain-outs, both scintillators detected a significant count enhancement, 
or a burst, over a 90 sec interval of 15:45:10 — 15:46:40 UT. Excluding the data obtained 
in this interval and applying a quadratic least square fitting to the remaining data, we 
estimated background (solid curves in Fig. [1]) in the two scintillators. Subtracting the in- 
terpolated background from the total observed counts in the 90 sec period and calculating 
the background errors using the derived quadratic function and the number of bins used in 
the fit procedure, 117, both Nal and plastic scintillators gave highly significant increases of 
1208 ± 150 (8.1(t) and 658 ± 104 (6.3(t), respectively. 

Figure [2] gives details of the count histories of the NaT scintillator, together with outputs 
of the electric field and light sensors. Since the plastic scintillator signals are used here 
in anticoincidence, the Nal events are thought to represent energetic photons. The excess 
counts in the 0.01 - 3 MeV and 3-12 MeV bands [Fig. EJ^a) and (b)] become 1033 ± 
159 and 91 ± 15, respectively, while those without anticoincidece in the two energy bands 
are 1030 ± 159 and 92 ± 16, respectively. Since the excess counts with and without the 
anticoincidence thus statistically agree, we infer that the Nal signals were dominated by 
photons, not charged particles. 

On this occasion, the indoor BGO detector also showed a moderate increase in > 2.6 
MeV energies, with a ~ 4o" significance when summed over a 15 sec interval. The absolute 
time of this increase is consistent with that in the outdoor detector, within a relatively large 
uncertainty of 3 min with which the clock of the data-taking computer of the indoor detector 
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was adjusted. Supposing that the outdoor and indoor detectors observed the same event, the 
lower significance in the latter, compared to the former, suggests that the burst source was 
located outside the building, and the radiation was attenuated by the roofs and/or walls. 

Due to a limited dynamic range of the electric field mill, the electric field [Fig. [2](c)] was 
saturated at —100 kVm^^ at 15:45:10 UT (burst onset), and kept negative polarity until 
15:46:29 UT. At that time, the electric field rapidly changed its polarity from negative to 
positive within 1 sec, and saturated at +100 kVm~^. However, the outputs of the optical 
sensor [Fig. [2](d)] remained at mV which is typical at midnight (15:00 UT corresponds to 
local midnight): its small fluctuations, ~ 5 mV, reflect the ADC resolution. In reality, when 
the optical sensor detects lightning, it typically records a few hundred mV or even > 5 V on 
rare occasions. Thus, we conclude that no lightning occurred during the burst, and hence 
this burst is associated with thunderclouds. 

Figure M^a) shows the background-subtracted photon spectrum, summed over the burst 
period (90 sec). Here, we subtracted a background spectrum averaged over 15:30 — 15:40 
UT and 15:50 — 16:00 UT, when the thunderstorm was ongoing. This is to remove effects 
of Radon rain-outs, which gradually increased the < 3 MeV background up to twice those 
in quiescent periods. The resultant background-subtracted spectrum exhibits a very hard 
continuum, clearly extending to 10 MeV. 

We can estimate the initial photon energy spectrum at the source, from the background- 
subtracted spectrum. Theoretically, the number of non-thermal bremsstrahlung photons 
per unit energy emitted by monochromatic electrons is in inverse proportion to the pho- 
ton energy, up to energies which are close to the electron energy 2l|]. Also, a numerical 
simulation 22] has shown that a population of rather fiat power-law distributed, or even 
monochromatic, relativistic electrons colliding with a thick medium produces a power-law 
bremsstrahlung photon spectrum with a very hard photon index. Thus, we simply assume 
that the initial photon spectrum is a power-law of the form aE~^, where a is a normal- 
ization factor (photons MeV^^sr^^), j3 is the photon index, and E is the photon energy in 
MeV. Then, to derive the photon energy spectra to be observed at the observatory, we sim- 



23j incorporating 



ulated the photon propagation in the atmosphere with CORSIKA 6.500 
EGS4 2J| to correctly treat electromagnetic interactions in relatively low energies relevant 



to the present work. The photons were assumed to be injected vertically into the atmo- 
sphere, and to propagate while suffering Compton scattering and other processes. Finally, 
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the simulated arrival photon spectra, after convolving with the detector response, were fit- 
ted to the observed photon spectrum in search for the best-fit values of a and /5 for various 
source distances assumed. The predicted spectra for representative source distances up to 
1000 m are shown in Fig. El^a); larger distances give poorer fits (e.g. = 42.2 for 3000 m). 
From Fig. [3]^b), the source distance, d = 90 m, gives the minimum /d.o.f. = 16.8/16, with 
a = (2.4 ± 0.4) X 10^ MeV~^sr"^ and (3 = 1.15 ± 0.09 (quoted errors are 90% confidence 
values). Figure [3]^b) also gives a constraint as = 60 — 130 m at 90% confidence level. 



DISCUSSION 



The present 7-ray burst is not due to Radon or its daughters, because the 7-rays related 
to Radon have energies of 0.1 — 2 MeV, and last for a much longer time tlO|. Compared with 

riQ 

some prolonged intracloud x-ray emissions [8|, |9|, this burst strongly suggests that electrons 
were accelerated beyond 10 MeV in thunderclouds. Furthermore, the most outstanding 

difference of the present burst from previous ones , 0, y, Q Q, Q, H, Q is that 

charged particles, most likely electrons, were simultaneously observed. This difference may 
be attributable to the short distance, ~ 90 m, of the present event. Actually, 20 (10) MeV 
electrons have a range of 110 (60) m at an altitude of 2770 m, which is in good agreement 
with the constraint on d derived above. 

We can construct the following picture of the present event. As predicted by Gurevich 
et al. [l7|, seed electrons with energies of > 100 keV, produced by e.g. cosmic rays, are elec- 



trostatically accelerated in thunderclouds to relativistic energies, and are multiplied therein. 
We may expect the electron spectrum to reach ~ 20 MeV, because a numerical study 25 1 
predicts that the runaway electrons have an average kinetic energy of ~ 7 MeV with a spread 
of ~ 12 MeV. The 20 MeV electrons can effectively emit bremsstrahlung 7 rays extending to 
10 MeV, and escape from the thunderclouds to propagate a distance of ~ 90 m, with their 
energies reducing down to a few MeV. Because of a Compton optical depth of 0.51 (0.12) at 
1 (10) MeV for d = 90 m, the emitted 7 rays are sometimes scattered by large angles with 
significant energy loss, producing secondary electrons. 

If the primary electrons were in a perfect parallel beam, relativistic effects would make 
the emitted 7 rays, especially high-energy ones, mostly beamed into a forward narrow cone 
with a half-opening ang le of F"^ ~ 1.5°, where F = 40 is the Lorentz factor of 20 MeV 
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electrons. In practice, the cone will be much broader, due to multiple scatterings of the 
emitting electrons, and to the mild Compton scatterings of the emitted photons. The 7- 
ray spectrum is expected to be hardest along the cone axis. Since the observed spectrum 
clearly extends to 10 MeV, it is thought that a rather harder part of the 7-ray emissions was 
detected. In particular, the derived photon index, 1.15 ± 0.09, is close to the theoretically 
hardest limit (1.0). Thus, we infer that our detectors viewed the cone nearly along its axis. 

We may estimate how the Compton-scattered electrons and the escaping primary ones 
contribute to the excess counts detected by the plastic scintillator, Noh = 658 ± 104. Ac- 
cording to a Monte Carlo simulation, the survival probability of > 1 MeV (at the detector) 
secondary electrons, produced by 3, 5, and 10 MeV 7 rays propagating from a source at 
d = 90 m, are 0.01, 0.02, and 0.04, respectively. Thus, employing the photon spectrum at 
d = 90 m, f9o{E) = 2.4 x lO^E'^-^^, the expected count of the > 1 MeV secondary electrons 
in the plastic scintillator becomes at most Uge ~ 110 even if all the scattered electrons with 
arriving energies of > 1 MeV are collected: we neglected lower-energy ones, since ioniza- 
tion losses in the 1 mm thick aluminum window prevents them from depositing > 500 keV 
energies in the plastic scintillator. Hence, we presume that the escaping primary electrons 
contribute to A^ob by Ad = Aob - '^sc = 548 ± 104. 

From this A^d, we can evaluate the total number of 20 MeV primary electrons for d = 90 
m as 

N20 = Ndd^iSe) - 3.1 X 10^ (1) 

where S = 1800 cm^ represents the area of the plastic scintillator, and e ~ 0.08 denotes the 
attenuation factor of > 1 MeV electrons in the atmosphere and the aluminum window. Using 
this A^2o and the initial photon spectrum fgQ{E), the spatial vertical length of the acceleration 
region, needed for the accelerated electrons to produce the observed bremsstrahlung 7 rays, 
can be estimated with an assumed beam half-opening angle 6h as 

H^27TN^o'j j f,oiE)/T]iE,e)smedEde, (2) 

where ri{E,6) is the probability per 1 gcm~^ with which the 20 MeV electrons emit 



bremsstrahlung photons with an energy E and an angle 6 with respect to the cone axis |21 |. 
At 9h = 15°, which is ten times larger than the perfect parallel beam case (1.5°), the above 
H becomes 200 m, and in turn gives the electrostatic potential difference in the accelera- 
tion region as U = 200 kVm-i X 200 m = 40 MV, where 200 kVm'^ is the critical electric 



field for seed electrons to cause the runaway electron avalanches given by 280P kVm~^ [20| 
with the atmospheric pressure P = 0.72 atm. This U is generally sufficient to accelerate 
high-energy seed electrons to relativistic energies or ~ 20 MeV if an electric field is higher 
than the critical value by 30 — 40%. As discussed so far, this observation suggests that 
long-duration emissions of 7 rays and electrons are due to relativistic runaway electrons, 
although other possibilities may not necessarily be excluded. 
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FIG. 1: Count rates per 30 seconds of > 10 keV Nal (top) and > 500 keV plastic (bottom) 
scintillators between 15:15 and 16:15 UT on 2008 September 20. The horizontal axis represents 
universal time, and all errors are statistical la. Each solid curve shows estimated background level. 
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FIG. 2: Detailed count histories of the NaT counter with anticoincidence, and the electric field and 
light sensors, between 15:40 and 15:52 UT. All abscissa are universal time. Panels (a) and (b) 

show count histories per 7 seconds of the Nal counter in 0.01 — 3 MeV and 3 — 12 MeV energies, 
respectively, with la statistical errors, (c) One-second electric field data variations, (d) One-second 
optical data outputs. 
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FIG. 3: (a) The background-subtracted photon energy spectra (black points) obtained from the 
Nal signals with anticoincidence. The horizontal and vertical axes represent the photon energy and 
the counts per energy bin, respectively. Error bars are statistical la. Each curve shows prediction 
of an incident power-law model. Air mass in g cm^^ corresponding to the assumed source distance 
is shown in parentheses, (b) The value of the model fit to the data, plotted against the assumed 
source distance. Horizontal lines from bottom to top represent 68%, 90%, and 99% confidence 
levels. 
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